Incommensurate Matrix Product State for Quantum Spin Systems 
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We introduce a matrix product state (MPS) with an incommensurate periodicity by applying 
the spin-rotation operator of each site to a uniform MPS in the thermodynamic limit. The spin 
rotations decrease the variational energy with accompanying translational symmetry breaking and 
the rotational symmetry breaking in the spin space even if the Hamiltonian has the both symmetries. 
The optimized pitch of rotational operator reflects the commensurate/incommensurate properties 
of spin-spin correlation functions in the 5* = 1/2 Heisenberg chain and the 5 = 1/2 ferro-antiferro 
zigzag chain. 
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I. INTRODUCTION 

An analysis of low-dimensional frustrated quantum 
spin systems beyond the mean-field approximation 
(MFA) is one of attractive topics in quantum mechanics, 
because rich quantum phases can appear due to the co- 
existence of frustration and strong quantum fluctuation. 
A typical example is the spin S = 1/2 ferro-antiferro (F- 
AF) zigzag Heisenberg/XXZ spin chain as a theoretical 
model of quasi-one dimensional edge-sharing cuprates. In 
theoretical studies on this quantum Hamiltonian^"^, the 
exact diagonalization method (ED), the density matrix 
renormalization group method (DMRG)^"— and the infi- 
nite time-evolving block decimation method (iTEBD)i^ 
were used as powerful methods in order to determine 
novel quantum phases. 

In the DMRG and the iTEBD methods, variational 
states take the form of a matrix product state (MPS)^'^.!^^ 
and an infinite MPS (iMPS)^, respectively. When the 
dimension of matrices constructing the MPS is one, the 
MPS corresponds to the MFA. As the dimension m in- 
creases, the optimum variational state approaches the ex- 
act one systematically. In addition, the MPS can handle 
infinite system-size directly if we suppose the spatial ho- 
mogeneity of the MPS as in the iTEBD. As a merit of the 
spatially uniform MPS or IMPS, there are no boundary 
effects which always appear in the DMRG. 

In the zigzag chaini"— , the helical magnetic order with 
incommensurate period is known to be a solution of the 
classical vector spin Heisenberg model which is valid in 
the large spin limit (S* 3> 1). The incommensurate prop- 
erties appear due to the geometrical frustration. To deal 
with quantum fiuctuation, one can use the MPS. How- 
ever, the spatially uniform MPS with finite dimension 
cannot express the helical magnetic order, because its 
local magnetic moment becomes spatially uniform. On 
the other hand, the DMRG can deal with a spatially in- 
homogeneous magnetic order, but the boundary affects 
incommensurate period of the order. 

In this study, we propose a simple incommensu- 
rate (IC) MPS with incommensurate periodicity ap- 
plying spin rotation operators^^ to the spatially uni- 



form MPS. This IC-MPS is understood naturally as 
a quantum generalization of the classical vector spin 
analysis. This framework is independent of the type 
of numerical optimization process, and it is applica- 
ble for various variational methods based on finite di- 
mensional MPSs: DMRG^^ii , the wave function pre- 
dictions based on the product wave function renormal- 
ization group (PWFRG) methodi^"— , the tensor prod- 
uct state (TPS )^ — ' ■— , the projected entangled pair state 
(PEPS)23, iTEBDi2, the infinite PEPS (iPEPS)Si, the 
tree tensor network (TTN) stated, the multiscale entan- 
glement renormalization ansatz (MERA) stated, and so 
on. To demonstrate our light-weight modification for the 
uniform MPS with small dimension of matrices m, the 
modified Powell method^! is used as a general purpose 
optimization method in this paper. 

A pitch angle which determines an incommensurate 
period is a variational parameter in our approach. The 
pitch angle plays an important role in the optimization 
of the variational energy. This is caused by the finite m 
effect, because any state can be expressed by the MPS 
with infinite m. However, in the analysis of quantum 
effect starting from the classical vector spin model our 
approach shows a fast convergence with respect to m 
and a result obtained by tiny m is consistent with IC 
spin-spin correlation propertiesii^. 

The spatial periodicity and translational symmetry are 
recent hot topics for the MPS and its generalizationSSr^. 
Our previous study'^^ shown that in the spatially uniform 
MPS the translational symmetry breaking appeared in 
principal eigenvalues of its transfer matrix; that is, the 
degeneracy of eigenvalues was consistent with the ground 
state periodicity. This means that we need large dimen- 
sion of matrices for the spatially uniform MPS with one- 
site periodicity to express a magnetic ordered state with 
p-site commensurate periodicity. To reduce computa- 
tional memory without loosing the numerical accuracy, 
p-site periodic MPS was effective^. However, as shown 
in this study, we succeed in reducing more computational 
memory using the IC-MPS. 

This paper is organized as follows. In ijlll we 
review the interaction-round-a-face (IRF) /vertex- type 
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]Yjpg32-34 g^^^ propose the IC-MPS. The MFA hmit of 
the IC-MPS is discussed in i^IIIl where we show that 
the optimum vertex-type IC-MPS with m = 1 in the 
5 = 1/2 Heisenberg chain is equivalent to the state from 
the MFA. In ^IVl observing the m dependence of local 
magnetization in the S* = 1/2 Heisenberg chain, we con- 
firm the IC-MPS takes into account the quantum fluc- 
tuation gradually by increasing m. In the same section, 
the effectiveness of the IC-MPS is demonstrated in the 
magnetization curve of 5' = 1/2 Heisenberg chain and 
the 5 = 1/2 F-AF zigzag chain under uniform magnetic 
fleld. Then, we discuss the reduced computational cost 
by applying the spin rotation in the 5 = 1/2 Heisenberg 
chain and the C-IC change with respect to the spin-spin 
correlation in the zigzag chair^""— i^ii^. 



II. MATRIX PRODUCT STATE WITH AN 
INCOMMENSURATE PERIOD 

Let us recall the IRF/vertex-type MPS^^^. An IRF- 
type MPS with TV site is 



Vertex-type IC-MPS 



(1) 



where Uj means the index of spin at ith site and a — 
ai ■ ■ ■ aj^. The variables Aq""^^ and ' are m x m 
square complex matrices. The matrix Aq"'^^ is called 
the boundary matrix^^iiliHi^. A vertex-type MPS is 
represented under the constraints: Aq"°^^ = A'^«Aq and 
A"^''^'^^ = A"^' . To handle the thermodynamic limit 
{N — !> C!o), hereafter, we treat a uniform MPS, namely 
ylj' = A'^i'^'+i. As in the previous study-—, one can 
treat the p-site periodic MPS. 

To construct an IC-MPS, we use a spin- rotational op- 
erator at each ith sitei^: 



Ri{m, Qi) = exp(-i(5iSj ■ rii), 



(2) 



where i means a unit of pure imaginary number and 
represents the local spin operator. The unit vector of 
rotational axis and angle at each site are represented by 
rii and Qi, respectively. In this paper, we limit ourselves 
to the simple case of = n and Qi = iQ. 
Then, the IC-MPS is given by 



I*) =i?tot(n,Q)|*). (3) 



A schematic picture of the wave function of IC-MPS is 
depicted in Fig. [l] where Ri means a matrix representa- 
tion of the operator R^. 

The variational energy for a Hamiltonian H is given 
by e(^, n, Q) = limAr^oo £'(*, n, Q)/iV with 

£;(*,n,Q) = (*,n,g|ifl*,n,Q)/(*,n,Q|^',n,Q} 

= {^\H{n,Qm/{^\^), (4) 
F(n,g) =i?L(n,Q)7?i?tot(n,g), (5) 
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FIG. 1. (Color online) Graphical representations of vertex- 
type IC-MPS and IRF-type IC-MPS. Filled circles and filled 
small squares mean contraction with respect to the local spin 
bases ai and the local artificial bases ai of the matrix A, 
respectively. 



where H[n, Q) is the spin-rotated Hamiltonian. Then, 
hereafter, we just consider H[n, Q). 

For general and Qi, the important characters of 
the spin-rotated operator are summarized below. The 
rotated local spin operator in general spin S is given by 

St{ni, Qi) = R\{ni,Qi)SiRi{ni, Qi) = T>{ni,Qi)si, (6) 

where the three dimensional matrix Ti(ni^Qi) is given 
by 



q)] = v^v^, + {S,^^, - v^v,y ) cos q 
-smq^e,„',y'V^„, 



(7) 



for the unit vector v. Symbols 3^^^^, and e^^'^/' represent 
the Kronecker delta and the Levi-Civita symbol, respec- 
tively, where rj = x,y, z. From Eq. ([7]), we can immedi- 
ately obtain the relation D(ni, Q^)' = D(ni, —Qi). 

For the simple case of = n and Qi = iQ, one can 
prove the following equation: 

Sj(n, iQ) ■ Si+e{n, [i + ()Q) = ■ Si+e{n, iQ). (8) 

The vanishing of position dependence simplifles the cal- 
culation of the Heisenberg Hamiltonian. For the 5=1/2 
Heisenberg chain defined by 



Si ■ s, 



+1' 



(9) 



the spin-rotated Hamiltonian Hi{n, Q) is written as 

Hi{n, Q) = Y,s,- s,:+i(n, Q) = ^ ^'("' Q)- (10) 

i i 

This Hamiltonian has the translational symmetry. Then, 
we apply the same uniform MPS used in the previous 
study—. If the artificial translational-symmetry break- 
ing does not occur, we can neglect the boundary matrix 
and the local energy = {^\hi{n,Q)\^) becomes inde- 
pendent of position i in the thermodynamic limit. The 
translational symmetry of the spin-rotated Hamiltonian 
is recovered even for the zigzag and bilinear-biquadratic 
Heisenberg chain for general spin 5. 
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It should be noted that we can deal with the case that 
the spin-rotated Hamiltonian does not have the trans- 
lational symmetry. In this case, we can calculate the 
variational energy by using the translational symmetry 
of the MPS 1 5'), because the position dependence of the 
local energy Ci can be expanded as 



= (0) 



3^=) exp(iifcg) 



(11) 



and only e^"-* gives non-zero contribution after taking the 
summation if Q is not commensurate. For com- 

mensurate Q, we must consider the contribution from 
e^*'-' for k ^ 0. Of course, one can treat more general 
position-dependent rotations, for example Q2i — '2-iQa 
and Q2i+i — (2i + l)Qbi where the expansion as in 
Eq. pT|) becomes more complex, namely = e'^"'"^ -I- 



(fca,fcb)#(0,0) 



=(fea,feb) 



exp(li(fca(5a + fcbQb))- 



0.01 
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FIG. 2. (Color online) Energy error AE as a function of local 
magnetization \M\ in the S — 1/2 Heisenberg chain Hi. The 
number nearby each symbol means the matrix dimension m. 
The broken line is guide for the eyes showing the power low 
decay of AE oc \M\'^. 



III. MEAN FIELD APPROXIMATION LIMIT 

We derive the mean-field limit of this method, which 
is realized by the vertex- type IC-MPS with m ~ 1. In 
this limit, we can neglect the boundary which has 
only trivial two roles: normalization and phase factor. 
Then, the MPS becomes a direct product state Y^) — 
Scr riil^'^' Wi))^ expressed by two complex variables, 
and va S = 1/2 systems. As a normalization, we 
assume Yla^ \A'^^\'^ = 1- 

To show that the mean-field limit corresponds to the 
classical vector spin model, we consider the Heisenberg 
Hamiltonian, Hi . The variational energy is given by 

e(*,n,g) = M-(D(n,Q)M), (12) 

with an expectation value of local magnetic moment 
M = I]^^/ ^'^*^'''(o'|s|c7'). The local magnetization is 

obtained by \M\ = After the optimiza- 

tion for fixed Q, one can obtain 

e(Q) = mine(5',n, Q) cosQ/4. (13) 

Then, the optimization of e{Q) gives the Neel-type so- 
lution Q = n. We stress again that this energy gain of 
e{Q) — e(0) is due to finite m, because any state can be 
expressed by the uniform {p = 1,Q — 0) MPS accurately 
if we have enough large dimension m for the MPS. This 
finite-dimensionality also causes \M\ = 1/2 which is al- 
ways proved for any state in the mean-field limit, while 
it is known that the exact ground state does not have 
the magnetization at zero magnetic field. In this sense, 
the mean-field limit corresponds to the classical vector 
spin model. In fact, as shown in mVl when we increase 
m to express quantum fiuctuation or entanglement, the 
local magnetization \M\ obtained after the optimization 
decreases and approaches to the exact value. 



IV. NUMERICAL RESULT AND DISCUSSION 

Before showing results, we summarize details of our nu- 
merical calculation. We prepare the m-dimensional com- 
plex matrix A'^"'^'+^ for the IRF-type uniform MPS. The 
rotational axis n is fixed as (0, 0, 1) to conserve transla- 
tional symmetry of the rotated uniaxial Hamiltonian with 
the longitudinal magnetic-field applied in z-axis. The 
pitch Q and A'^^-'^'+i are optimized so that the variational 
energy e for a given Hamiltonian H becomes minimum 
by using the modified Powell method^-. The number of 
optimization parameters in the IRF-type IC-MPS under 
fixed rotational axis is 2(Pm? + 1, where the coefficient 2 
comes from using complex numbers and d is the degree 
of freedom (DOF) of local spin, namely 2 in this work. 
The term g-l-lh means the DOF of the wave number Q. 
In the optimization, 10 — 2000 initial states are prepared 
and optimized in each Hamiltonian parameter to avoid 
obtaining a local minimum. 

The MPS gradually takes account of the quantum 
fiuctuation of the local magnetic moment in the S=l/2 
Heisenberg chain with increasing m as shown in Fig. [2] 
The rotational angle Q — it is, obtained after the op- 
timization. The energy error /S.E means the difference 
between optimized variational energies as function of m 
and the exact energy — In 2 + 1 /4'^^ . The energy error and 
the local magnetization \M\ are monotonically decreas- 
ing with respect to m. We confirm that the IRF-type 
MPS can deal with non-zero quantum fiuctuations even 
if TO = 1, while the vertex- type MPS with to = 1 gives 
the mean-field result. This is an advantage of using IRF- 
type MPS. 

The magnetization Mz in the = 1/2 Heisenberg 
chain with the magnetic-field is shown in Fig. [31 As 
reference data, we show the exact result for S — 1/2 from 
the Bethe ansata2£ and the result from two-site modu- 
lated MPS, named p = 2, Q = 0, in our previous study—. 
While the mean field result fails to obtain the correct crit- 
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FIG. 3. (Color online) Magnetization curve as a function FIG. 4. (Color online) Rotational parameter Q dependence 
of uniform magnetic field in the 5 = 1/2 Heisenberg chain of the variational energy in the S = 1/2 Heisenberg chain Hi. 
Hi. 



icality near the fully saturated point, results for m = 3, 
which is not so large dimension, show enough accuracy. 

This increasing of m leads to a great improvement in 
the accuracy of estimating the magnetization curve. The 
relative error of the magnetization curve from the IC- 
MPS with m = 3 in Fig. [3] is smaller than 3% even 
though the error of local magnetization \M\ is of the 
order of that from the MFA as shown in Fig. [51 that is, 
the absolute error of Mz in Fig. ^is less than 0.001 even 
though the error of \M\ in Fig. [2] is about 0.1. 

Moreover, in both cases of m = 1 and m — 3, the 
data from IC-MPS with Q = tt agree with that of the 
p = 2,Q — MPS. This means that the number of opti- 
mization parameter is reduced by 50% compared to the 
previous study. 

How the rotational pitch Q = tt is stabilized by the 
energy gain, e{Q) — e(0), is shown in Fig. |4l This fig- 
ure clearly shows that the variational energy becomes 
minimum at Q = tt for any magnetic field except for 
Hz ~ 2.0 in the perfect-ferro region. Compared with 
= cos (5/4 for = in the mean-field limit, there 
is the fiat energy region in small Q region for — 0. 
In the flat region, we conflrm the state is a superposi- 
tion of the Neel state, namely the linear combination of 
1 tit • ■ •) and I ttl • • •)^^- This state is invariant with 
respect to the spin rotation along z-axis trivially. The 
origin of the flat region is the quantum fluctuation of the 
Neel state. This quantum fluctuation can be expressed 
by the IRF even in m = 1. 

Finally, we discuss the periodicity change appearing in 
the S = 1/2 F-AF ( Ji < and J2 > 0) zigzag Heisenberg 
chain with uniform longitudinal magnetic held, 

i \k=\.:i I 

The longitudinal magnetic field Hz is taken as and 0.1 
in this analysis. At Hz — 0, there is a C-IC change 
at JxjJi = —4. The commensurate state for J1/J2 < 
—4 is a ferromagnetic state while characterization of the 



ground state for J1/J2 > —4 is a difficult task. Recent 
study^ pointed out that the ground state for J1/J2 > 
—4 is the Haldane-dimer phase, which is characterized 
by a generalized string order parameter, where ordinal 
spin-spin correlations behave incommensurately'^. This 
incommensurate behavior is also found in the VC phase 
for non-zero magnetic fieldai. 

To demonstrate our approach for the C-IC change, op- 
timized pitch Q is calculated for this frustrated Hamil- 
tonian Hi as shown in Fig. [5] In this figure there 
are three kinds of the reference data. First, the bro- 
ken line is the result of the mean-field approximation, 
Q = arccos(— J1/4J2). Second, the solid line is the fitting 
line for the location of the maximum of the zero field spin 
structure factor with the ED3., where Q oc (J2- l/4)°-2^. 
Finally, the filled circles are the result of the DMRG at 
Mz = 0.05 in the VC phased For Hz = 0, the pitch Q 
approaches 7r/2 with increasing J1/J2 more rapidly than 
that of the mean field approximation due to taking into 
account the quantum fluctuation by growing m. On the 
other hand, the C-IC change point is completely con- 
verged at J1/J2 = —4. We find the pitch in m = 3 is 
comparable with the result from ED"^ in J1/J2 < —2.5. 
Around the transition point, the pitch is well converged 
with respect to m in this scale. The m dependence be- 
comes gradually large with increasing J1/J2, where the 
frustration due to J2 becomes also gradually large. 

For Hz = 0.1, the pitch Q depicted by the open circle 
in Fig. [5] has a jump around Ji/ J2 = —3.1, which is very 
close the SDW3-VC phase transition pointi. Unfortu- 
nately, the pitch Q for small m fails to capture the SDW3 
and SDW2 states, where there are the Ferro-SDWa phase 
transition at J1/J2 ~ —3.3 and the VC-SDW2 phase 
transition at J1/J2 —2^. In the same meaning, the 
characterization of the ground state at Hz = is dif- 
ficult for our method at this stage. Nevertheless, a no- 
table point is that the incommensurate pitch Q of the IC- 
MPS for the VC phase in J1/J2 > —3.1 shows reasonable 
agreement with the DMRG result. For J1/J2 > —3.1, Q 
of the IC-MPS is nearly independent of Hz, which is also 
consistent with the DMRG analysis^. 



5 




FIG. 5. (Color online) Rotational parameter Q as a function 
of J1/J2 in the zigzag chain i/2. 



V. SUMMARY 

In summary, we introduced the IRF-type MPS with 
the incommensurate pitch parameter Q and the rota- 
tional axis n as a generalization of the uniform MPS, 
which can be used for various variational methods based 
on the MPS. Two parameters Q and n allow us to evalu- 
ate an incommensurability of the spin chain in the ther- 
modynamic limit directly. Our approach with small di- 
mension of matrices is connected to the classical vector 
spin Heisenberg model which is valid in the large spin 
limit [S ^ 1). For the exact ground state, the helical 
magnetic order obtained in the classical limit is expected 
to be destroyed by quantum fluctuations in the quantum 
limit S = 1/2. However, we emphasize quantum effects 
on some quantities are rapidly converged with respect 
to the matrix dimension. Our approach opens a way 
to a light-weight analysis based on the classical vector 
spin model to include quantum fluctuation. Using this 
approach, one can treat translational symmetry broken 
states, such as the helical magnetic order, in the thermo- 
dynamic limit, which cannot be handled by known IMPS 
with translational symmetry. 

We demonstrated the efficiency of this IRF-type IC- 
MPS in two types of Hamiltonians: i) the magnetization 
in the S = 1/2 antiferro- magnetic Heisenberg chain un- 
der uniform magnetic field, and ii) the C-IC change in 
S = 1/2 F-AF Heisenberg zigzag chain under uniform 
magnetic field. In the former Hamiltonian, we have suc- 
ceeded in obtaining the same result as two-site modulated 
MPS. This means 50% reduction of the number of opti- 
mization parameters. In the latter Hamiltonian, we have 
succeeded in detection of the C-IC change of correlation 
properties with increasing to. The pitch Q near the C-IC 
transition point is immediately converged with respect to 
TO and shows reasonable agreement with the ED study'^ 
and the DMRG studjsi, despite small to. 

On the other hand, the sufficiently converged Q is not 
obtained around the strongly frustrated region, namely 
|Ji| ~ J2. To discuss the details of Q, analysis with 
larger m are necessary. For this problem, we can ap- 



ply other optimization methods using the Trotter de- 
composition-i^, the matrix product operator (MPO)^, 
and the time-dependent variational principle (TDVP)^ 
to updating the MPS under given Qi and n^. We stress 
again that the framework of IC-MPS is independent of 
the type of numerical optimization process. In this pa- 
per, the modified Powell method was chosen as an opti- 
mization method because it is a general purpose method 
and all parameters are optimized easily. The modified 
Powell method is enough to clarify to effectiveness of our 
light-weight modification but becomes a bottleneck when 
we increase to. To study larger m, convergence proper- 
ties and numerical efficiencies of these updating methods 
should be discussed. This is one of future problems. 

Another future issue is to change the constraint of the 
rotational axis and pitch parameter in order to represent 
the magnetization plateau state or the SDW state, for 
example rii = n^jod^j and Qi — Qniod[i.p]- As another 
application, we have already performed other C-IC corre- 
lation properties change in the bilinear-biquadratic spin 
S — 1 chain", and succeeded to detect the C-IC change 
with the IC-MPS, which cannot be detected by the mean 
filed approximation. 

This method uses the spin-rotational operator which 
maps the classical helical state to the perfect-ferro state. 
The uniform direct product state including the perfect- 
ferro state can be always described by the uniform MPS 
with to = 1. A generalization of this method is to find 
another kind of spin-rotational operator which maps the 
ground state to the uniform direct product state. The 
role of this operation is similar to disentanglers in the 
MERA^S. In this sense, it is interesting to consider the 
valence bond solid (VBS) state which cannot be rotated 
by the spin-rotational operators^- and is known to have 
the Kennedy- Tasaki (KT) transformation which converts 
the string order to the ferromagnetic order as a global 
topological disentangler^^'^-. 

In general, the classical magnetic order can be ap- 
peared easily in higher- dimensional systems. In this case, 
the spin rotation becomes effective. Moreover, in higher- 
dimensional systems, the dimension of the matrix/tensor 
is restricted due to the computational resources. Then, a 
small TO analysis based on our approach is an interesting 
approach for the incommensurate TPS for 2D quantum 
spin systems, which is one of future issues. Not only for 
the TPS, our approach can be applied to various meth- 
ods. 
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